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CYCLES TILL FAILURE OF SILVER-ZINC CELLS WITH COMPETING FAILURE 

MODES - PRELIMINARY DATA ANALYSIS 

by Steven M. Sldlk, Harold F. Lelbcck'l^ and John M. Bozek 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


INTRODUCTION 

The sllver/zlnc electrochemical couple has found many applica- 
tions. Mainly, Its use In the past has been limited to primary ap- 
plications requiring high currents, low weight, and permitted 
limited charge/discharge capabilities. With the advent of space 
probes and orbital spacecraft, secondary applications also became 
Important. These applications require high specific energy and 
charge /discharge cycle life but permit limited current capability. 

A sealed secondary sllver/zlnc cell was developed for this reason 
(ref. 1). This cell used an Inorganlc/organlc (I/O) separator and 
was very successful In the laboratory (refs. 2, 3, and 4). It was 
rated at 40 ampere hours (Ah), delivered 88 w-hr/kg (40 w-hr/pound) 
and about 400 cycles at a 40 percent depth of discharge (refs. 3 
and 4). 

Some near earth orbit missions motivated the development of a 
12 Ah version of this cell. The specifications of this cell are: 

12 Ah rated sealed cell 

Pressed silver powder positive electrodes 

Pressed zinc oxide powder negative electrodes 

I/O Bag Separator around both electrodes 

Three negative electrodes 

Four positive electrodes 

45 weight percent KOH electrolyte 

Size: 

12.0 cm high to top of case 
6.0 cm wide 
2.25 cm thick 
316 gm weight 

The specific energy of this cell was judged superior to that of 
other candidate electro-chemical couples (such as Nl/Cd, Ag/Cd) 
since the Ag/Zn couple is the most active. The theoretical spe- 
cific energy of the silver-zinc cell is 440 w-hrs/kg, of the nickel- 
cadmium cell is 210 w-hrs/kg and of the silver-cadmium cell Is 270 
w-hrs/kg. It remained, therefore, to characterize the charge/dis- 
charge cycle life of this silver-zinc cell. 
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An experiment wee deelgned which would serve two purposes. The prim- 
ary objective of the experiment was to characterize the cycle life per- 
formance as a function of five variables (use conditions). These five 
Variables are 


Charge rate (CR) 

Discharge rate (DR) 

Depth of discharge (DOD) 

Ambient temperature (T) 

End of charge voltage (ECV) 

A second objective of the experiment is to estimate the variability of cycle 
life among cells tested ui$der nominally the same conditions. 

This report will first present a discussion of the variables considered 
and the experimental design. The design is a variation of a central cempos- 
Ite factorial design commonly used In response surface methodology (ref. 5). 

We then turn tb fitting cycle life as a second order polynomial func- 
tion of the five variables and conclusions that can be drawn from this fit. 

It will be shown that the primary results are the strong effect of depth 
of discharge (DOD) on cycle life and that cycle life peaks at about 25 ° to 
30° C ambient temperature. 

In the following section we differentiate between the two primary fail- 
ure modes* shorting and low voltage. We first describe a competing failure 
modes model and estimation procedure (refs. 6 and 7). We then present equa- 
tions which predict cycle life until a short develops if a low voltage cannot 
occur and cycle life until a low voltage failure If a short cannot occur. 
Whether a battery first falls by low voltage or shorting Is determined pri- 
marily by temperature. At the lower temperatures the cells fall by low 
voltage and yield low cycle lives. At the highest temperatures the cells 
fall by shorting and yield low cycle lives. At Intermediate temperatures 
we observe a mix of failure modes and higher cycle lives. 


EXPERIMENT DESCRIPTION 


Experiment Design 


For this program, a total of 129 cells were tested. All were construc- 
ted from nominally Identical materials by nominally Identical methods. The 
five Independent variables and the ranges over which they were Investigated 
are 


Charge rate (CR) 

Discharge rats (DR) 

Depth of Discharge (DCH3) 
Temperature (T) 

End of Charge Voltage (E CV) 


0.375 to 1.625 amperes 
1.25 to 5.00 amperes 
20% to 100% of rated Ah 
0° C to 40° C 

1.98 volts per cell to 2.02 volts 
per cell 
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The particular combinations of levels of these variables were chosen as a 
variation on central composite factorial designs (ref. 5). Table I pro- 
vides the data for each cell tested. The test conditions are given in 
columns 2-6, the cycles to failure in columns 12-15, and mode of failure 
In column 16. 

In this experiment one power supply unit was used to cycle a number 
of cells collected in series at the same charge rate and the same dis- 
charge rate. Nine combinations of charge and discharge rate were inves- 
tigated. These levels are given in Table II. Each combination required 
a separate power supply unit. 

Within each combination of these levels a full 2^ factorial with 3 
center point replicates was run with respect to D(X), T, and ECV. The 
specific levels are given in Table III. In addition to these runs, the 
axial points of a 2^ central composite design on DOD, T, and ECV were run 
on power supply units 2, 4, 5, 6, and 8. The specific levels are also 
given in Table III. 

These tests were performed by the Naval Weapons Support Center, 
Crane, Indiana. 


Formation and Cycling Procedure 

Before subjecting the cells to cycling under their respective test 
conditions, each cell was put through a pre-test formation and capacity 
check. This Involved two charge/discharge cycles at room ambient temper- 
ature. Each cell was charged at 0.27 amps to 1.98-2.00 volts or 13.5 A-H 
whichever occurred first. Then each cell was discharged at 1.13 amps to 
1.3 volts followed by a drain at 0.377 amps to 1.3 volts. This charge/ 
discharge sequence was followed twice. The ampere-hours In and out of 
the cell on each of these charges and discharges were recorded. For 
purposes of later analysis, the capacity (A-H) removed from the cells at 
the 1.13 amp and 0.377 amp rates were added and simply called capacity. 
Also for purposes of later analysis, the capacity removed for both forma- 
tion cycles was averaged and Is reported In Table I as actual capacity 
(A-Hgctual) • 

After this formation the cells were ready for life cycling tests. 

Life cycling consists of discharging at the stated DR until the Ampere- 
hours out of the cell equals (DODjjom) (12)/100. If the cell cannot supply 
this at a voltage of 1.3 volts per cell or more, we say a low voltage 
failure occurs at that point. After discharge, the cells are charged at 
the stated CR until the voltage across the cells reaches ECV. If this 
takes more than twice the time It should at the stated charge rate we say 
a failure occurs. From the charge/discharge voltage behavior and the 
ability of a cell to hold a charge It Is evident when a cell has developed 
an internal short. A cell may, therefore, fall by one of two modes, low 
voltage or shorting. 
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The complete cycling till failure process la as follows. After the 
formation» cells are cycled until a failure condition occurs. If this 
failure condition is a short, testing terminates and the cycle at which 
this occurs is denoted f (in Table I) . If this first failure condition 
is a low voltage condition, the cycle is noted as fj_ in Table I for 
that cell but charging and discharging continue as planned until the 
secovnd failure condition occurs. If this second failure condition is a 
short the cycle at which it occurs is denoted in Table I and further 
testing is terminated. If this second failure condition is another low 
voltage condition, the cycle at which it occurs is denoted £.2 in Table 1 . 
The cell is then reformed once according to the formation criteria of 
voltage and current outlined previously. 

After this reforming, cycling is continued until a third failure con-* 
dition is encountered. If this is a short the cycle is recorded as f^ 
in Table I, and testing terminates, wlaile if it is a low voltage condition 
we note the cycle as 13 in Table I and continue cycling as planned. The 
fourth failure condition is denoted as f^ in Table I and terminates life 
cycling. 

Since shorting terminates testing whenever it occurs, there are no 
further cycle lives reported in Table I after a short occurs. Refonuing 
of the cells takes place only after cycle f2- 


ANALYSIS OF OVERALL CYCLE LIFE 

In this section we consider life (number of cycles till failure) as 
a response regardless of the failure mode. In the following section we 
differentiate between the two modes of failure and describe equations for 
predicting the number of cycles until a cell fails by one mode under the 
condition the other mode is Inoperative. 

The five independent variables we use for life prediction are; 

» charge rate (CR in amps) 

^2 “ discharge rate (DR in amps) 

^3 ® depth of discharge (DOD in X) 

5^ " ambient temperature (T in °C) 

55 = end of charge voltage (ECV in volts) 

In terms of these variables, characterising cell performance is equivalent 
to providing an equation predicting cycles till faiL-;re as a function of 
53, ..., 55. Rather than do this directly we define scaled variables (as 
is common in response surface methodology (ref. 5 )) for the obse:-va« 

tion as 
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hi “ 

^2j ” ^hi " 3 *13) A -87 

hi " “ 67.2)/19.4 (1) 

hi • (^4j “ 20) /lO 

hi “ ‘ 2.00)/0.OX 

For X,, Xnf Xa, Xe the center point value Is subtracted and the increment 
from the center point to the high level is used as the divisor. For X 3 , 
the nominal depth of discharge (DODj^on,) as indicated in Table I is not 
used but the actual depth as calculated from the actual capacity based 
upon the formation cycles is. That is, 


C 


3 


12 DOD 


nom 


^ " *^actual 


( 2 ) 


For X 3 , the center point of 67.2 is the average of 5 - while the scaling 
factor of 19.4 is the standard deviation of 63 , 

The nominal DOD is the value a designer might use for design or missioa 
analysis purposes where the actual capacity and actual DOD are generally un- 
known. But, once the actual cell is on hand and put through a capacity 
check, its actual capacity and hence actual DOD may be determined for pre- 
dictive purposes. 

As the response variable we use Yj = Yj(Xj^, ..., X 5 ) which may be cycles 
till failure (f^^, 62 , or f^) or some function of cycles till failure 
(e.g., log fi). 

The candidate function chosen as a starting point of analysis is simply a 
second order polynomial as commonly used in response surface methodology. 

That is, the statistical model is 


Vj . + 8iX„ + • • • + 65X5^ 

+ Pixxf^ 


(3) 


+ 615X1JX5J + ... + + e 

2 

where e is a random error with zero mean and constant variance 0 and the 
3 's are unknown coefficients to be evaluated empirically. 
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As the response variable we chose to use y ■ logj^Q (cycles ) . This 
choice was first made by analysis of residuals from fitting using y " cycles 
and supported by probability plots of the residuals using y * logj^Q<cycles) . 

For each cell there are up to four failure cycles (depending upon when 
a short occurs). For purpose of analysis , if a cell failed by shorting before 
cycle f/, we simply repeat the last observed failure cycle for all succeeding 
failures? For example, cell number 626 failed first by the low voltage mode 
at cycle 74 and then failed by shOi ting at cycle 75 (see Table 1) . For a 
shorted cell, reconditioning is not possible and we simply use 75 for f, 
and f^. There would be a qualitative difference between a cell which shorted 
at cycle 75 before ever failing by low voltage and say a cell which failed by 
low voltage twice before finally falling by shorting at cycle 75. This would 
not cause any serious distortion in modeling f3 failures since the fact still 
remains that both cells failed by the 75th cycle. 

Two cycle life predictive equations were fit, one to predict the £2 fail- 
ure cycle and the other to predict the f^ failure cycle. The rationale for 
this is as follows. When the cell has a low voltage condition at cycle f^ 
that indicates an inability to meet mission specifications for that cycle. Yet, 
it is possible to obtain more cycles which meet specs after that and Indeed, 
Table I shows this to happen in almost every low voltage situation. Essenti- 
ally, this is because low voltage failure is a slow degradation phenomenon and 
it is difficult to pinpoint a precise point of failure. Calling £2 the life 
of a cell is a somewhat arbitrary attempt at defining such a point of failure. 
The cells are then reconditioned and again tested till "failure" which we will 
consider to be cycle f^. The first equation, then attempts to predict useful 
cycle life if reforming cannot be done while the second equation attempts to 
predict useful cycle life if up to one reformation is permitted. 

Cells number 618 and 727 were not Included for analysis since they failed 
due to different causes (i.e. , cell 618 failed due to operator error and 727 
failed when its case ruptured). 

Variable X5 (measuring the effect of ECV) had no significant main effect 
or interaction with any variable in any of the analyses performed and is not 
Included in any further discussion or analysis. (This lack of significance is 
probably due to the restricted range over which ECV was varied) . 

Table IV provides the summary statistics and estimated coefficients of 
equation (3) for the £2 cycle life. The first column Identifies the coeffi- 
cient of equation (3) whose estimate is being reported. The second column 
gives the estimated coefficients (estimated by linear least squares regression 
analysis) and the estimated standard errors of the coefficients enclosed in 
parentheses. Those coefficient estimates significantly different from zero at 
the five percent significance level are indicated by an asterisk. The stand- 
ard errors and significance level are determined by pooling both the pure 
replication and the lack of fit errors. At the bottom of the table are given 
the number of data points used in the fitting (n = 127) ^ the standard error 
of estimate (S = 0.342) and the multiple correlation coefficient squared 
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(r2 m 0,789). The third column of Table V gives the estimates and standard 
errors of the coefficients when the data is refitted using only those coef- 
ficients significant at the five percent level of significance (i.e., re- 
duced model). The standard error of estimate (S « 0.348) and squared 
correlation coefficient (r2 w 0.771) indicate that the two f equations fit 
the data essentially Identically and will provide very nearly equal pre- 
dictions. 

Both model equations indicate that there are only two significant inter- 
actions among the variables (the interaction between charge rate and depth 
of discharge and the interaction between discharge rate and temperature). 

The Interaction between discharge rate and temperature means that the partic- 
ular effect temperature has on the cycle life of a cell is dependent upon 
the rate at which it is discharged. 

Table V gives the corresponding information for the f^ failure. Com- 
parison of Tables IV and V shows very little difference either qualitatively 
or quantitatively between the fitted equations. The only appreciable differ- 
ence is a lower standard error of estimate reported in Table V. 

These equations and data were examined for Indications of lack-of-fit or 
spurious observations by examining various residual and probability plots. On 
the basis of these it seemed that the results of cells 722 and 726 did not fit 
the patterns. Similar probability and residual plots discussed in the nex't 
section also indicated that cells 602 and 608 were not consistent with the 
rest of the data. Thus, equation (3) was refit to the 123 cells remaining 
after deleting from the analysis cells 602, 608, 722, and 726. The fitted 
constants for equation (3) and summary statistics for f2 failures are given 
in Table VI and the results for failures in Table VII. Comparison of 

the full equations using 127 cells as opposed to 123 cells show t'ae equations 
to have quite similar coefficient estimates. When comparing the reduced 
versions of the equations (i.e., all coefficients significant at the five 
percent level) there are minor differences in the second order terms. Spec- 
ifically, for f2 failures (Tables IV and VI) the X^, X§, and X^X^ terms 
are different. For failures (Tables V and VII) the XjlX 4 interaction 

appears in the 123 cell fit but not in the 127 cell fit. For the terms which 
appear in both full equations the coefficients are very similar in magnitude. 
The major changes are in the degree of fit as evidenced by the smaller stand- 
ard error of estimate and larger r 2 values for the edited data set (i.e., 
the data set containing only 123 points) . 

We examine the fitted response surfaces in two ways. First, we perform 
a canonical analysis of the full second order equation. Second, we provide 
several parametric plots of life versus the independent variables. 

In a canonical analysis of a second order polynomial response surface we 
first determine the stationary point of the surface and second perform a rota- 
tion of axes about that stationary point to new independent variables. In this 
process we let IC^ denote the stationary point, Yg denote the value of Y 
at that point, and denote the new (rotated) axes. (For details and further 
descriptions of this type analysis see ref. 5.) This reduces equation (3) to 
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the form of equation (A) 


V -V3 - AizJ+ + XpZj 


(« 


where the are the eigenvalues of the matrix of second order coeffic- 
ients. The variables are the axes of symmetry and are given by 


^i " ®il^^l " ^Is^ + ’’* e^^pCXp - Xpg) 


(5) 


where the e. . are components of the eignvectors of the matrix of second 
order coefficients. 

Table VIII provides the parameters of Interest in the canonical analy- 
sis for the £2 cycle life and Table IX for the £^ cycle life. In 
Table VIII the Indicate that the response surface is a saddle-shaped 
one. The value X^^ » 0.057 Indicates a rising ridge type of structure 
with respect to the variable 

Zi = - 0.21(Xi + 4.67) + 0.18(X2 + 2.50) 

+ 0.95(X3 - 0.50) + 0.16(X4 + 1.17) 

The primary component of this variable is X 3 (depth of discharge) . The 
value X 4 “ - 0.274 indicates a falling valley with respect to the 
variable 

Z 4 « - O.nCX^ + 4.67) - 0.49(X2 + 2.50) 

- 0.09(X3 - 0.50) + 0.85CX4 + 1.17) 

This new variable is primarily a combination of X 2 (discharge rate) and 
X^ (temperature) . (This arises from the X 2 X 4 interaction which is large 
and significant as seen in Table VI.) The surface is almost constant with 
respect to variable Z 2 and somewhat decreasing with variable Z3. Veri- 
able Z3 is given by 

Z 3 » - 0.69(Xj^ + 4.67) + 0.62(X2 + 2.50) 

- 0.30(X3 - 0.50) + 0.19(X4 + 1.17) 

which is a combination of all four original variables. 

The X^ and e^^ of Table IX are very similar to those of Table VIII. 

The stationary point is shifted more toward the center of the design. Thus, 

the basic shape of the response surface is the same but shifted in location 
slightly. 

The shape of the surface may be easier to understand by looking at fig- 
ures 1 through 5. In these figures we plot the predicted logiQ^Z cycle 

life using the reduced equation of Table VI. Figure 1 plots *2 versus POD 
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for five different temperature values. For these plots CR is held constant 
at 1.0 amp and DR is held constant at 3.13 amp. This plot clearly shows the 
rising ridge structure with respect to temperature and DOD. Figure 2 plots 
f 2 versus D(^ for three charge rates. The small Interaction between CR 
and DOD shows a steeper loss of life with respect to DOD at the higher 
charge rates. The curves show that slower charging gives better cycle life. 
Figure 3 shows f 2 versus DOD at three discharge rates. There is no lnter~ 
action between these variables and the conclusion to be drawn is chat lower 
discharge rates give longer life. Figure A plots f 2 versus temperature 
at three different discharge rates. We see that maximal life is obtained 
at a temperature dependent on discharge rate. This implies that as the dis- 
charge rate required increases* the optimal temperature also rises. Fig- 
ure 5 plots f 2 versus temperature at three different charge rates. The 
Interaction between temperature and charge rate is not as pronounced as the 
interaction between temperature and discharge rate. For any given CR there 
is a temperature providing maximal life and this temperature Increases as 
charge rate Increases. The less pronounced interaction of figure 5 as opposed 
to figure 4 is shown by the three curves being more nearly parallel. 


competing failure modes analysis 

In the preceedlng section we developed equations for predicting failure 
cycle as a function of the test conditions (i.e., the independent variables). 
In this section we differentiate between the two failure modes low voltage 
and shorting. Before the cell is put on test we assume a random time until 
failure due to low voltage has a probability distribution F(l)(y]^; X) and 
Independently a random time until failure due to shorting has a probability 
distribution F( 2 )(y 2 ; X). The actual life of the cell is the smaller of 
these two lifetimes. As the notation indicates, we allow the distributions 
F(i)(y£; X) to depend upon the test conditons, X. Specifically, we assume 
that F(i) ■= F(i)(yjLJ p(i»j), a (i)) where F is an extreme value distribu- 
Vion with location parameter and scale parameter The location 

parameter is assumed to be related to the test conditions (similarly to eq. 
(3)) as 

p(ij) = g(i) + + ••• + 


+ S^^^x2 

11 Ij 


+ X 4- 

®12 + hz 




(i) 2 
+ B 55 


( 6 ) 
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For the Bmallest extreme value distributions 

F(y; MfO) " 1 - exp exp | (7) 

(If life follows a Wcibull distribution, then log life Is known to follow 
a smallest extreme value distribution.) The observed data is of the form 


life mode 


test conditions — - 



“l Hi ’ ' ‘ ^Ik 

« • • 

i f • 



M 


n 



« f « 



We also critically examine the assumptions using probability plots and 
"residual” plots. 


Methods for estimating the parameters (i.e., the B’s and o's) of 
these models and some simulation studies of the properties of these estim- 
ators is given in reference 7. The method used is maximum likelihood and 
this requires the solution of sets of non-linear equations. (The equations 
to be solved are given in the appendix.) The method of solution is straight- 
forward but can be unstable if the model equation (eq. (6)) has too many 
terms which are insignificant. In an attempt to avoid this potential prob- 
lem, we first separate the full set of data into two parts. The first part 
of the data is for cells that failed by low voltage and the second part for 
cells that shorted. Models of the form of equation (3) are then fit using 
ordinary least squares to each data set individually. This provides a 
first guess as to the form of the equation needed for each mode and also 
starting values for the g*s and o's required for the iterative maximum 
likelihood method. 


As with ordinary multiple linear regression analysis, these models must 
be considered tentative in nature and examined for adequacy of fit and 
spurious observations. We present our results by discussing: 


- The least squares fits to the short-only and low voltage-only 
data 


- The corresponding maximum likelihood estimates 

- Various "residual" plots for diagnosing adequacy of fit 

- And the maximum likelihood estimates for the models and edited 
data sets 
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Analyses of All the Dato 

The data was split into sets of * 97 low voltage failures and 
n« ■ 30 short failures for a total of 127 dato points. As discussed 
earlier » two cells (727 and 618) have been eliminated from the analysis. 

Table X presents the results of the linear regression analysis for cycle 
£2 failures while Table XI presents the corresponding results for cycle 
f 4 failures. In each table we retain only those terms with coefficients 
significant of the five percent significance level except for the T and 
t 2 terms. For the short data there is a large correlation between T 
and t 2 when using the scaling of equations (1). It turns out that either 
one but not both of these terms are significant. We Included both of those 
terras regardless of their individual significance levels. Tables X and XI 
do not require much discussion except to note that the values of the g's 
and or's reported therein were used as starting values for the iterative 
maximum likelihood method. 

In Table XII we present the maximum likelihood estimates of the coef- 
ficients for the f 2 failure cycle for each mode and their estimated stand- 
ard errors in parentheses. Upon comparing Tables X and XII it may be seen 
that the coefficient estimates are not far apart except for a few instances. 

For the low voltage failures > the DR coefficients are somewhat different. 

For the short failures the individual coefficients of T and T^ are some- 
what different. For both modes the constant coefficients are larger by the 
maximum likelihood method than by the least squares method. Similar comments 
hold true for the f^ failure results (reported in Table XIII) except that 
the XxX2 interaction has reversed signs for the short failure mode between 
Tables XI and XIII. 

These fits were examined by "residual" plots as follows. For the separ- 
ate least squares fits, standardized residuals were calculated in the standard 
manner (l.e., (y - y/s). These were plotted in various ways including versus 
each Independent variable and versus, the predicted value. Three of these plots 
are given in figure 6 through 8. Figure 6 presents the residuals against dis- 
charge rate for the cycle f 2 low voltage failure regression. There is some 
Indication that the variance increases as DR increases and that there may be 
three outlying cells at the highest DR. (These are indicated by the O and 2 
at the bottom right of the plot.) Two of these three are cells 722 and 726. 
Figure 7 presents residuals against scaled temperature for the cycle f^ low 
voltage failure data. The two largest negative residuals are cells 722 and 
726. There is also a large positive residual at the lowest temperature. Fig- 
ure 8 plots the residuals versus the predicted values for the cycle low 

voltage failure regression. It almost appears as if there are two populations 
of residuals. The lowest predicted values are over-predicted while the resid- 
uals tend toward more uniform behavior as y increases. This probably indi- 
cates some sort of flaw in the model equation. One problem with the model as 
used here is that for the vetry low cycle lives, logiQ (cycles) cannot be well 
described by any continuous distribution, such as the extreme-value or normal. 
The greatest number of cells had much larger lifetimes, however, and this is 
where we are most interested in battery performance anyhow. The only action 
taken at this point was to delete cells 722 and 726 from analysis. 
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From the equations estimated by maximum likelihood and presented in 
Tables XII and XXIX ve computed "residuals" as follows. We first compute 
p(x) the predicted location parameter for the mode by which that cell 
failed. The residual is then defined as y^^ - u. This is an expedient 
but not necessarily correct approach and may lead to some distortion where 
there is considerable overlap of the distributions of the two failure modes. 
Va.rious plots of these residuals were examined and three of these are given 
as figures 9 through 11. Figure 9 plots the residuals versus scaled charge 
rats for the cycle f^ short failures. The two large negative residuals 
at the lowest charge rate are cells 602 and 608. 

In figure 10 we provide a probability plot of the residuals on an 
extreme value probability scale while figure 11 gives a probability plot 
On a normal probability scale. These are for the cycle f^ short failures. 
It is quite clear that the residuals fit an extreme value distribution much 
better than a normal distribution. The only action taken at this point was 
to delete the data from cells 602 and 608 from the analysis. 


Analysis of Edited Data 

The data were re-analysed by the same techniques described in the pre- 
vious sections except with the data from cells 602, 608, 722, and 726 removed. 
Equation (3) was fit to the 93 low voltage failures and the 28 short failures 
separately with the resulting least squares coefficients given in Tables XIV 
and XV. Comparing Table X which is for the unedited data to Table XIV 

which gives the fit for the edited f 2 data, it may be noted that there is 
very little difference between the coefficients. The fits seem somewhat 
Improved as evidenced by the larger r 2 and smaller S values. Similar 
comments hold for the comparison of Tables XI and XV which present the f^ 
regression results, 

As before, the least squares results were used as initial values for the 
iterative maximum likelihood procedure. Table XVI gives the refined estimates 
for the l>i failure times. In parentheses next to each coefficient estimate 
are given their estimated •standard errors. Comparing Tables XII and XVI we 
see that the charges in coefficient estimates are minor except for the indiv- 
idual coefficients of T and T^ for the shorting mode equations. Comparing 
Tables XIII and XVII which give the maximum likelihood estimates for the 64 
cycle life, the differences are again all small ejtcept for the Individual 
coefficients on T and T^ and the reversal in sign of the DR x CR inter- 
action. The sign and magnitude of the DR x CR interaction in Table XVII are 
consistent with the least squares results. 

It may be noted that in Tables XVI and XVII there are several coefficient 
estimates which are within one estimated standard error of zero (e.g., the £2 
DR low voltage and f 2 DOD^ short coefficients of Table XVI) . One is tempted 
to compute a t-like statistic and conclude these terms are not significant 
This may not be true because correlations among the 0 and I estimates may 
lead to a very skewed dlstribuLion of the t-ratio in some Instances (ref. 8 ). 
We have refitted the maximum likelihood estimates for reduced models but these 
are not presented or discussed in this report. It turns out that the coeffi- 
cients of Tables XVI and XVII provide the lowest (most conservative) estiraatori 
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for life. For this reason we use these equations to compare the effects of 
the independent variables on the life by each modeo Figures 12 through 16 
provide some of the more Interesting plots of predicted life times. 

For a smallest extreme value distribution with location parameter p 
and scale parameter o, the expected life is p - where y “ 0.5776 
is Eulers constant.^ For the following plots, we use as expected life 
p(Xi» ••• , X 5 ) - ya. 

Figure 12 plots expected (log life £ 2 ) versus temperature for three 
different discharge rates. For these, charge rate is held at 1.0 amp, 

DOD at 67.2 percent of actual capacity and E CV limit at 2.00 volts. The 
curves for low voltage failure cycle indicate initially rising life with 
temperature. There is a maximum with respect to temperature and the max- 
imizing temperature depends upon discharge rate. The higher the temperature 
and discharge rate, the more cycles before a low voltage failure occurs. The 
curves for short life indicate that as temperature increases, shorts develop 
sooner and that shorts develop sooner at the low discharge rates. At about 
25° to 30° C is the crossover between failure modes. At the lowest tempera- 
tures the cells fail by low voltage long before shorts have had a chance to 

develop. At 40° C, the cells are developing shorts before they have the 

chance to fail by low voltage. One implication of the reversed behavior 
of life with respect to temperature is that if a change could be made to 
the cell which prevented low voltage failures from occurring at tempera- 
tures of say 10° C but had no effect on shorts developing, then we could 
expect an Increased life of approximately 400 cycles. It is not clear that 
such changes could be made. Certainly, more effort should be expended on 
post-mortem an/.lyses of cells to break down cause of failure into more 
directly assignable causes than low voltage and short. 

Figure 13 plots expected (logj^Q life £ 2 ) versus depth of discharge 
for several different temperatures; For these, charge rate was held at 

1.0 amp, discharge rate was held at 3.13 amp, and EOCV limit at 2.00 volts. 

The plots show that expected life decreases with depth of discharge but at 
slightly different rates. The dependence of failure mode on temperature 
is clearly evident also. 

Figure 14 plots logiA (life f^) versus DOD for various charge rates. 
For these, temperature is ■‘'held at 20° C, discharge rate at 3.13 amps, and 
E CV limit at 2.00 volts. For the low voltage failure mode, expected life 
is decreasing with DOD and, at these conditions, is the dominant failure 
mode regardless of charge rate. For the shorting failure mode, expected 
life decreases with DOD at the higher charge rates while it essentially 
remains constant at the lower charge rates. 

Figure 15 plots log-.^ (life f^) versus DOD for various discharge 
rates. For these, charge^rate is held at 1.0 amp, temperature at 20° C 
and E CV limit at 2.00 volts. For the low voltage failures, there is 
little effect on life from discharge rate but expected life is decreasing 
with DOD. For the shorting mode expected l° 8 io (f 2 ) decreases with DOD 
at the low discharge rates but decreases less at the higher discharge rates. 
As in figure 14, the dominating failure rode is low voltage at those condi- 
tions . 


Figure 16 plots expected cycles till S 2 Failure for three charge 
rates. It Is similar In shape to figure 12 which plots expected cycles 
till f 2 failure for three discharge rates. Figure 16 shows that at low 
temperatures low voltage Is the dominating failure mode while at high 
temperatures shorting Is the dominating failure mode. 


DISCUSSION 

In this section we discuss the effects of each of the independent 
variables, the results of the regression on overall life and their inter- 
pretation, the results of the competing failure mode modal and analysis, 
and recommendations for further research. 

There are numerous factors which could affect the cycling life of a 
cell but the five we considered are primarily those that define the use 
environment of the cell. Their levels and ranges of variation were chosen 
on the basis of what was anticipated to effect cycle life at the time the 
experiment was designed. 


Charge and Discharge Rates 

Variations In charge and discharge rates have an effect on life for 
several reasons. Either rate may have an effect on the crystal structure 
of the active materials which may in turn affect the rates at which the 
electrochemical couple progresses. High charge rates may have accelerated 
zinc dendrite growth leading to early shorting (ref. 9). In fact, fig- 
ures 2 and 5 show that overall life decreases as the charge rate Increases. 
Figure 3 shows that at T = 20° C the overall life decreases as DR Increases. 
Figure 4 shows that at low temperatures overall life decreases with DR while 
it Increases with DR at higher temperatures. Figure 12 indicates the effects 
by mode. That Is, at the lowest temperatures, the dominating failure mode 
Is low voltage and the effect of DR on low voltage life Is mixed. At the 
higher temperatures, the dominating failure mode is shorting and shorts 
develop In fewer cycles at the lower discharge rates. 

Figure 5 shows that overall life decreases as CR increases. Figure 16 
shows this to be true for each mode« It is interesting that decreasing the 
discharge rate promotes early shorting while increasing charge rate promotes 
early shorting. Thus, the data Indicates slow charging and rapid discharg- 
ing are best in terms of shorting. 


Depth of Discharge 

Variations in the depth of discharge require different percentages of 
the active materials to be exercised during each cycle. Since the reactions 
are never totally reversible, high DC© should lead to decreased cycle life, 
In fact, figure 1 shows a large decrease in overall cycle life. Figures 13, 
14, and 15 also show that both low voltage and shorting lifetimes generally 
decrease as DOD Increases. At the lowest charge rates and the highest dis- 
charge rates there is little or no decrease in cycle life as DOD increases. 
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Care must bo taken in drawing conclusions at coiulltlons other than specified 
in the figures due to the internctiono among the variables. 


Temperature 

Rates of chemical reactions depend strongly upon temper.-Jturo and hence 
cycle life of a coll should depend strongly on the temperature. Changes in 
temperature may effect the crystal growth, charge acceptance, growth of 
dendrites, or dissolution of cell separator components. Figures 4 and 5 
show the effect of temperature on overall cycle life. For each charge and 
discharge rate combination there is some optimal temperature. Figures 12 
and 16 sJiow the mode of failure depends strongly on the temperature. High 
temperatures promote j'vrocosses that lead to shorting while at low tempera- 
tures shorts are slow to develop. 

The rate of zinc dendrite propagation Increases as the concentration 
of the zlncate ions increases while the initiation time for dendrite form- 
ation decreases with increasing over-potentJnl , temperature and zlneate 
concentration (ref. 9). At Xo\i temperature, the diffusion of zlncate is 
less than at the higher temperatures and would cause a mean transport con- 
ditions conducive to the results observed. 


End of Charge Voltage 

It was felt that E CV might effect cycle life primarily due to the 
possibility of gassing at high end of charge voltages. In this experiment 
we found no evidence of an effect of E CV over the range examined. It is 
possible tlut we did not investigate a largo enough range. 


Life Prediction Models 

For mission analysis and component design purposes, the analyst is 
primarily interested in knowing what cycle life might be expected under 
the conditions at which the cell might be used. In general, the mode of 
failure is unimportant. (This is to some extent not true because failure 
by low voltage is not a well-defined condition. For most missions a short 
failure may be considered to be terminal whereas if low voltage failures 
occur, the mission might bo extended if the cell is reformed.) The equa- 
tions of Tables VI and Vll provide predictive equations for overall 
failure by our definitions of failure. The reliability of the equations 
as a predictive tool are partially indicated by the estimated standard 
deviations which imply that the uncertainty In a further predicted 
logfo (f2) is S ■ 0.300, 


Plus and minus two standard error prediction limits in loRin ^^2^ 

2S “ 0.60 translate to error limits of fjj x I 2o or from 0.^5 f 2 to 

3.98 failure equation we havt S » 0.248 which leads to 

proportional error limits of from 0.32 to 3.13 f^^. That is, we can 
predict cycle life to wlcMn a factor of 3 or 4. 


For the battery researcher, the competing failure mode model and 
analysis provide considerable more Information for Improving cell cycle 
life. In this study we have simplified the category of failure Into two 
classes of failure, namely low voltage and shorting. This Is recognized 
as a simplification since there may be several places and causes for 
either type of failure. Ideally, each cell should undergo a post-mortem 
analysis to further differentiate falliue causes. The competing failure 
modes analysis does provide a first attvn.pt at separating causes and an 
Indication of which design changes might be most Important to consider. 

For instance, figure 12 shows that at a cemperature of 0*^ C and DR of 
1.0 amp a cell will yield about four cycles before a low voltage fail- 
ure. If a design change could be utilized which alleviated this type of 
failure but would not affect the formation of shorts, we could expect to 
get about 500 cycles before shorting would terminate cell life. Of 
course, such changes may also adversely affect cycle life due to shorting. 
The Independence of the two failure modes assumed In the model could be 
a restrictive assumption. However, the use of the competing failure 
mode model analysis along with post-mortem analysis of the cells 
promises to be a useful and informative tool. 
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APPENDIX 

MAXIMUM LIKELIHOOD ESTIMATE EQUATIONS 

Assume that an experiment is performed in which I items are life 
tested at varying combinations of J stress variables and that the design 
is specified by the design matrix 


^11 


'Id 


Lzii • • ‘ 


(A-l) 


The response observed for each item is a lifetime y^_ (i « 1*1) and the mode 
by which the item failed, mi* The number of modes by which the item might 
fall is denoted M. 


We assume that for observation i and mode m the cumulative distrib- 
ution function of time until failure is defined by 


where 




1 - exp 



(A-2) 


P 


(m) 

i 




'IJ 


(A-3) 


The probability density function is given by 



(A-4) 


In general we will not necessarily have 


gCm) B 0(m’) 


nor 


(A-5) 
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p(m) , g(m') fQy ^ 

3 3 

That io, we do not a priori assume any of the failure modes to have param- 
eter values In common. 

It may be shown that the natural logarithm of the likelihood function 
Is 



To obtain maximum likelihood estimators, the standard method (assuming it 
works) is to solve the system of JM (nonlinear) equations 


3 In L 
36 « 




« 0 (m = 1, M) 


3 In L 


83 


(m) 


= 0 (m = 1, M) 




(A-9) 


3 In L 


30 


(m) 


= 0 (m = 1, M) 


If the different failure modes have no parameters in common, as will gen- 
erally be true by equations (A-5) and (A-6) , this system of JM equations 
splits into M separate systems of J equations which may be independently 
solved. In particular, the m^*^ system reduces to 
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From equations (A-3) and (A-4) it may be shown these equations reduce 


to 
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1 

^ij “a 


failed by 
mode 
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These nonlinear equations have no closed-form solution but may be solved by 
several Iterative methods. 

It Is well-known that, under appropriate regularity conditions, maKlrnuro 
likelihood estimators (MLE) are consistent and asymptotically normally dis- 
tributed. It is also known for the Welbull and extreme value distributions 
that the minimal sufficient statistics are the trivial ones consisting of 
the order statistics. It thus becomes of considerable interest to determine 
how rapidly the MLE's approach their asymptotic unbiasedness and covariance 
structure. The asjrmptotic covariance structure and the Cramer-Rao lower 
bound are defined by the Inverse of the Fisher information matrix (ref. 10, 
p. 194). 

For the model considered in this report the Fisher information matrix is 
the block-diagonal matrix 


where 


l(l) 0 
0 I<2) 


l(M) 


(A-12) 
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1 
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(A-13) 


These expressions are estimated as described in reference 6 and yield the 
estimated standard errors of the coefficients reported in Tables XII, XIII, 
XVI, and XVII. 
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trABLE 12 AKP-HE SltVER-ZINC CHARACTERIZATIOW TEST 


PACK 

NO. 

AOIA 
402A 
A05A 
60AA 
A05A 
A06A 
607A 
(OBA 
60?A 
610A 
6UA 
B12A 
B13A 
BIAA 
615A 
6I6A 
6 17 A 
6UA 
619A 
A20A 
“*621A 
622A 
623A 
62AA 
"Z625A 


CHARGE DISCHARGE DEPTH OF 
RATE RATE DISCHARGE 

AMPS AMPS X 

TE51 

TEST 

* CONDZTK 
END OF 
CHAKOE 
V01T5 

ACTUAL FORMATION 
CNOIl DCHGOl CHGI2 
A-H A-H A-H 

CAPACITY 
DCH0I2 CH003 
A-H A-H 

— — TEI 
Cycles to 

PI F2 

5T RESULTS 
FAILURE 
F3 F4 

FAILURE 

MODE 

0,375 

1*25 

03.2 

10 

i.ft 

11. 89 

10.55 

11.76 

11.22 

11.66 

138 

195 

157 

U7 

IV 

0.375 

1.25 

70. B 

10 

l.ff 

11.74 

10.62 

11*39 

11.02 

11.43 

101 

103 

121 

122 

s 

0*375 

1.25 

03*2 

30 

1.90 

12.61 

11*71 

11,71 

11.13 

12.27 

296 

297 



s 

0.375 

1*25 

70. B 

30 

1.99 

I2.lt 

10, *4 

U,87 

11.30 

11*99 

150 

151 



s 

0.375 

1.25 

03.2 

10 

2.01 

12.29 

11.55 

12,25 

11.60 

12.42 

187 

201 

290 

291 

IV 

0*375 

1*25 

70. B 

10 

2.01 

12.45 

11.21 

11.26 

10.85 

12*02 

76 

7^ 

92 

93 

LV 

0.375 

1.25 

03.2 

30 

2,01 

12. B7 

12.11 

12.31 

U.44 

12.78 

295 

252 



s 

0.575 

1.25 

70. B 

30 

2*01 

13*18 

11,96 

n.79 

10,97 

12,05 

86 




s 

0*375 

1.25 

00.0 

20 

2.00 

11.99 

10.69 

11 .91 

10,96 

12*09 

325 

331 

337 

338 

LV 

0.375 

1,25 

00.0 

20 

2,00 

11.70 

10*91 

11.97 

11.06 

12.02 

281 

282 

306 

308 

LV 

0.375 

1.25 

00.0 

20 

2,00 

13.48 

12.20 

12.60 

11,39 

12.81 

309 

310 

359 

360 

LV 

1*000 

1.25 

03.2 

10 

1.99 

12.99 

11.73 

U.61 

10.90 

12,92 

78 

79 

81 


IV 

i.ooo 

1.25 

70.B 

10 

1*99 

12*99 

11.70 

11.63 

10.38 

12.00 

8 




LV 

1,000 

1.25 

03.2 

30 

1.99 

12.83 

11.73 

11.67 

a. 09 

12,52 

283 

285 

292 

299 

LV 

1,000 

1.25 

70.8 

30 

2.01 

12.63 

11*60 

12,09 

a, 97 

12.30 

130 

131 



S 

1.000 

1,25 

03,2 

10 

2.01 

12.29 

11.93 

11.15 

10,59 

a. 88 

79 




LV 

1.000 

1.25 

71.3 

10 

2.01 

12.38 

11.36 

12,08 

a. 59 

12*15 

9 

92 

59 

55 

LV 

1.000 

1.25 

03.2 

30 

2.01 

12.29 

11.99 

11.86 

10,76 

12.96 

115 




OP 

1.000 

1.25 

70.8 

30 

2.00 

12.83 

U.69 

12.07 

11.99 

12.69 

160 

161 

169 

165 

LV 

1.000 

1.25 

00.0 

20 

2,00 

11*95 

10.26 

a. 12 

10.32 

a. 65 

122 

123 

UO 

169 

IV 

1.000 

1.25 

00.0 

20 

2.00 

11.17 

9.98 

a. 27 

10.79 

11*66 

193 

199 

159 

155 

LV 

1.000 

1.25 

00.0 

20 

2.00 

12*19 

11*11 

12.15 

a *87 

12.53 

m 

162 

172 

173 

LV 

1.000 

1.25 

20.0 

20 

2.00 

12.69 

11.87 

a. 23 

10,53 

12,95 

966 

967 



LV 

1.000 

1.25 

100.0 

20 

2.00 

12.20 

10,77 

10.98 

10*57 

11,66 

99 

100 

108 

110 

LV 

1.000 

1.25 

60.0 

0 

2,00 

11.07 

10.17 

a . 02 

10.95 

11,98 

3 

9 

5 

8 

LV 


TABLE I. - ContinMd. 12 AKP'-lffi 8ILVER-ZINC CHARACTERIZATION TEST 


TEST RESULTS 


fhtK 

HOi 

CHARCE 

RATE 

AflP5 

DISCHARGE 

RATE 

AMPS 

DEPTH OF 
DISCHARGE 
X 

TEST 

'fr 

END OF 
CHARGE 
VOLTS 

ACTUAL FORMATION 
CHGOl DCHOOl CMGI2 
A-H A-H A-H 

CAPACITY 
DCHOI2 CHG83 
A‘M A»H 

CYCLES TO 

FJ n 

FAILURE 
F3 F. 

FAUU) 

MODE 

42IA 

1.000 

1,25 

00.0 

00 

2,00 

10,81 

10,51 

11,20 

10,87 

11,15 

70 

n 



5 

I27A 

r.ooo 

1.25 

00,9 

20 

1,98 

11,01 

10.30 

11,00 

11,13 

11,75 

105 

106 

130 

135 

LV 

I2IA 

1.000 

1,25 

00.0 

20 

2.02 

11.55 

10,23 

11.00 

11.07 

12,02 

150 

177 

183 

180 

LV 

iZik 

1.125 

1.25 

03,2 

10 

1,99 

12,90 

11.83 

11,80 

10,89 

12,01 

55 

56 

75 

76 

LV 

i30A 

1.025 

1.25 

70,8 

10 

1,99 

12,05 

10,89 

11.70 

11.00 

12.51 

3 

0 

7 

8 

LV 

43U 

1.025 

1.25 

03.2 

30 

1.99 

11,97 

11.05 

10,80 

10.00 

11,08 

210 

217 



LV 

A 32 A 

1.025 

1.25 

70.8 

30 

1,99 

11,88 

10.00 

11,03 

10.28 

11,73 

02 

63 

80 

81 

LV 

4S3A 

1.025 

1,25 

03.2 

10 

2.01 

12,29 

11,17 

10.90 

10,13 

11.88 

79 

81 

96 

97 

LV 

A34A 

1,025 

1.25 

70.8 

10 

2.01 

11.85 

10,70 

11.00 

10,87 

12,10 

0 

5 

8 

9 

LV 

I33A 

1,025 

1.25 

03,2 

30 

2,01 

12,29 

11.38 

10.90 

10.17 

11.90 

229 

230 

238 

202 

LV 

I3AA 

1.025 

1.25 

70,8 

30 

2,01 

11.81 

10,00 

11.72 

10.98 

11,90 

109 

110 

113 

UO 

LV 

437A 

1.025 

1.25 

00.0 

20 

2,00 

12,02 

11,03 

11,03 

10.00 

12,20 

75 

70 

98 

99 

IV 

iSAA 

1.025 

1,25 

00.0 

20 

2,00 

11,01 

10,05 

11.21 

10,75 

12,21 

09 

70 

95 

90 

LV 

434A 

1.025 

1.25 

00,0 

20 

2.00 

11,97 

10.08 

11.90 

11,22 

12,20 

32 

33 

83 

80 

LV 

440A 

0.375 

3.13 

03.2 

10 

1,99 

12,29 

11.07 

11.35 

10.07 

12.00 

102 

170 

210 

200 

LV 

AAIA 

0.375 

3,13 

70.8 

10 

1.99 

12.30 

11.17 

11.03 

10.98 

12.02 

00 

50 

66 

71 

LV 

A42A 

0.375 

3.13 

03.2 

30 

1.99 

12,91 

11.92 

11,20 

10,53 

12.05 

398 

399 



S 

643A 

0.375 

3.13 

70.8 

30 

1,99 

11,23 

10,00 

10.99 

10.00 

11.15 

310 

315 

318 

319 

LV 

6A4A 

0.375 

3.13 

03.2 

10 

2.01 

n.oi 

10.61 

11 .48 

10.88 

12,31 

108 

109 

172 

170 

LV 

645A 

0.375 

3.13 

70,8 

10 

2.01 

12,10 

10.55 

11,33 

10.92 

11.95 

52 

53 

70 

71 

LV 

'646A 

0.375 

3. 13 

03.2 

30 

2.01 

12.00 

11.22 

11.08 

10.85 

12,11 

286 

287 



S 

A47A 

0.375 

3.13 

70.8 

30 

2.01 

11,03 

10.08 

11,53 

11,39 

11,01 

376 

377 



5 

64SA 

0.375 

3.13 

00.0 

20 

2.00 

12,08 

10.79 

11.00 

10.98 

12,15 

179 

180 

210 

219 

LV 

A49A 

0.375 

3.13 

00.0 

20 

2.00 

11.91 

10.02 

11.08 

b.l9 

11.16 

170 

171 

219 

221 

LV 

:a5oa 

0.375 

3.13 

00.0 

20 

2.00 

10.62 

9.00 

10.60 

9,83 

10,75 

205 

250 

275 

281 

LV 
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PACK 

HO, 

HMHMHMIlMNHNIIMKIItNtilMiHMIlNMItNMM TEST 
CHARGE DISCHAA06 DEPTH OF TEST 
RATE RATE 0I9CHAR0E TETIP 

AMPS AMPS X *C 

COMDITIOHI 
END OF 
CHARGE 
V0LT5 

ACTUAL PORMATION 
CHOU OCHOIl CH0I2 
A-H A-H A-H 

..K.K. ...... 

CAPACITY 
DCHGI2 CHQI9 
A“H A-H 

— — — TEI 

cycles to 
FI F2 

5T REIULtS 
FAILURE 
F3 F6 

FAILURE 

MODE 

ink 

0d375 

3.13 

20.0 

20 

2.00 

13.20 

11.06 

U.97 

10.3A 

11.85 

527 




5 

(S2A 

0.379 

3.13 

100.0 

20 

2.00 

11.61 

10,32 

11.10 

10,36 

11,08 

66 

65 

77 

71 

IV 

k$U 

0.379 

3,13 

60.0 

0 

2.00 

11,33 

10.06 

10.17 

10,00 

11.69 

1 

2 

3 

6 

IV 

454A 

0,379 

li.13 

60.0 

60 

2.0 

11.01 

10.66 

U.JA 

10.66 

11.79 

101 

102 



5 


0,375 

3.13 

60,0 

20 

1.06 

11,61 

10.23 

11.31 

11.06 

11.72 

107 

198 

217 

218 

LV 

kHk 

0,375 

3.13 

60.0 

20 

2,02 

U.,2 

10.23 

10,13 

0,«3 

10,20 

200 

203 

216 

217 

LV 

$57A 

1. 000 

3.13 

63.2 

10 

1.00 

11.77 

10.68 

U.«l 

10.00 

11.03 

33 

38 

56 

60 

LV 

i3H 

1.000 

3,13 

76.6 

10 

1,00 

U.IS 

11.13 

11,38 

10,66 

11,07 

16 

20 

26 

26 

IV 

an 

1.000 

3.13 

63,2 

30 

l.OO 

12.02 

10,60 

10.67 

0,08 

11.20 

217 

218 



5 

itOA 

1,000 

3.13 

76.1 

SO 

1,00 

U.77 

10.28 

11.66 

11,00 

11,03 

161 

162 



S 

Ulk 

1.000 

3.13 

63.2 

10 

2.01 

11.36 

0.08 

U,16 

10.62 

11.33 

6 

23 

25 

27 

LV 

HZk 

1,000 

3.13 

76,6 

10 

2.01 

11.36 

10.02 

11.31 

10.87 

11.99 

12 

13 

26 

27 

IV 

ki^k 

1.000 

3,13 

63.2 

30 

2.01 

13.32 

12.0, 

12.62 

11.26 

12.55 

360 




5 

Hkk 

1,000 

3.13 

76.6 

30 

2.01 

11.07 

10.06 

10,63 

10.26 

10.76 

187 

168 



5 


U000 

3.13 

60,0 

20 

2.00 

11,56 

10,63 


10,87 

12.16 

80 

100 

130 

U1 

IV 

mk 

1.000 

3.13 

60.0 

20 

2.00 

12.06 

10. M 

11.08 

10.28 

11,66 

06 

too 

210 

228 

LV 

kiik 

1.000 

3,13 

60.0 

20 

2.00 

11.50 

10.32 

11.35 

11.17 

U.56 

256 

262 

282 

283 

LV 

6 68 A 

1.000 

3.13 

20.0 

20 

2.00 

11.82 

10.55 

11.51 

10.68 

11.62 

089 

990 



S 

HH 

1.000 

3,13 

100. 0 

20 

2.00 

11.23 

10.26 

10.06 

10.70 

11.09 

68 

50 

67 

68 

LV 

670A 

1,000 

3.13 

60.0 

0 

2.00 

U 32 

10.17 

11.35 

10.68 

11.37 

2 

3 

6 

5 

LV 

“•671 A 

l.OOO 

3.13 

60,0 

60 

2.00 

U.IO 

10.11 

11.27 

10.75 

11.27 

127 

126 



S 

;‘H72A 

t.ooo 

3.13 

60.0 

20 

1.06 

11.17 

10.00 

11.16 

10.76 

11.09 

125 

126 

166 

169 

LV 

' '675A 

1.000 

3.13 

60,0 

20 

2.02 

11,56 

10.06 

10,72 

IP. 63 

11.33 

69 

75 

103 

106 

LV 

67AA 

1.625 

3.13 

63.2 

10 

1.00 

11 . 33 

10.17 

11,66 

10.51 

11.87 

35 

36 

53 

56 

LV 

*675A 

1.625 

3,13 

76.6 

10 

1.00 

11.63 

10.17 

10.62 

0.57 

U.87 

2 

3 

5 

6 

LV 


TABLE I. - Continued. 12 AKP-HR SILVBR-ZINC CUARACTERlZATIOK TE8T 


PACK 

NPe 

CHAA08 

f^ATE 

AHP5 

TEST 

pischarob depth op test 

RATE DISCHAROe TEMP 

AMP5 X *e 

’ CONDITIONS 

END OF ACIDAt rORMATlON 

CHAROE CHOIl DCHOIt CNa,2 

vons A-H A-H A-H 

CAPACirr 
DCKdl2 CH0I5 
A-H A-K 

TES 

CYCLES TO 
FI P2 

n RESULTS 
FAILURE 
F3 F4 

FAILURE 

HOPE 


i*m 

3.13 

43.2 

30 

1,99 

15.34 

12.07 

12.58 

U,21 

12.24 

305 




S 

♦7 7 A 

1.423 

3.13 

76, S 

50 

1.99 

12.14 

10,51 

U44 

11.09 

12,22 

159 

160 



S 

f7IA 

1.A25 

3.13 

43.2 

10 

2,01 

11,18 

10,45 

11.40 

10,51 

12,21 

43 

44 

48 

70 

lv 

*7»A 

1.423 

3.13 

76.1 

10 

2,01 

U«U 

10,02 

10,98 

1041 

10.78 

2 

3 

6 

7 

LV 

4B0A 

1.425 

3 43 

43.2 

30 

2,01 

11,32 

10.11 

11.09 

10.24 

1146 

279 

280 



$ 

ASIA 

Itl23 

3.13 

76. a 

30 

2,«t 

10.02 

9,15 

11.20 

10.32 

11,26 

148 

149 



5 

AS2A 

1.125 

3.13 

60,0 

20 

2.00 

11,63 

10. S5 

11,36 

10.62 

1146 

15 

16 

26 

106 

LV 

6aSA 

1.A25 

3.13 

60.0 

20 

2.00 

11,04 

9.64 

10.87 

10,00 

10.79 

13 

14 

18 

19 

IV 

iS4A 

1.425 

3.13 

60. 0 

20 

2.00 

11.84 

10.55 

11.91 

10.28 

11.99 

16 

19 

22 

62 

LV 

AS$A 

1.425 

3.13 

20.0 

20 

2.00 

10.89 

9.85 

11.19 

10,40 

10.94 

335 

357 

507 

508 

LV 

A86A 

1.125 

3.15 

100.0 

20 

2.00 

11,80 

10.68 

10.82 

10.24 

10,89 

27 

28 

40 

43 

LV 

A87A 

1.425 

3.15 

40.0 

0 

2.00 

11.69 

10.36 

11.56 

10.70 

11.75 

1 

2 

3 

4 

LV 

fiSSA 

1.625 

5.13 

40.0 

40 

2.00 

11.03 

1047 

11.51 

10,92 

11.56 

90 

91 



5 

6S9A 

1.625 

343 

60,0 

20 

1,98 

13.50 

12,13 

12.56 

11.30 

12,56 

43 

125 

135 


LV 

690A 

1.625 

3.13 

60,0 

20 

2,02 

12.41 

U.i5 

11.17 

10,49 

12.35 

188 

189 

232 

280 

LV 

691A 

54375 

5.00 

43.2 

10 

1.99 

12.72 

U.66 

11.80 

10,94 

12,56 

12 

15 

21 

22 

LV 

A92A 

0.575 

5.00 

76.8 

10 

1.99 

12.55 

U.4I 

11,48 

11.09 

12,55 

10 

12 

19 

20 

LV 

6 93 A 

0,575 

5.00 

43.2 

30 

1.99 

12,22 

11.25 

11.71 

10.74 

12.28 

503 




5 

A94A 

0.375 

5.00 

76.8 

30 

1.99 

11.51 

1047 

11.58 

10.40 

11,87 

159 

161 

163 

164 

LV 

A95A 

0*375 

9.00 

43,2 

10 

2.01 

12.0? 

11.04 

11,50 

10.26 

12.04 

15 

16 

22 

23 

LV 

'696A 

0.575 

5.00 

76.8 

10 

2.01 

12.02 

10.62 

11.40 

10.87 

12.40 

10 

12 

20 

21 

LV 

697A 

0.375 

5.00 

43.2 

50 

2.01 

12.55 

11.55 

11.71 

10.66 

11.94 

471 

472 



5 

69SA 

0.575 

5.00 

76,8 

30 

2.01 

11.82 

10,58 

11.86 

11.07 

11.75 

160 

161 

166 

167 

LV 

6 99 A 

0.375 

5.00 

60.0 

20 

2.00 

12.25 

u.n 

11, 4Z 

1041 

11 ,93 

103 

105 

134 

137 

LV 

'700A 

0.375 

5.00 

60.0 

20 

2.00 

11.20 

10.26 

11 V.8 

11.06 

11.97 

145 

147 

172 

175 

LV 


r 


'1 


TABU !♦ Concluded. 12 AMP-HK SILVrRR-74K0 CHARAOtKR! NATION imT 


PACK 

NO. 

TEST CONDITI0H5 

CHARGE DISCHARGE DEPTH OF TEST EHD OF ACTOAl FORMATION CAPACITY 

?ilf discharge temp charge CHUt) DCHOM CHOIE 0CH0*8 CHOU 

AMPS AMPS R »C V01T5 A-H A-H A"H A-H A-H 

CYCLES TO FAILURE 
FI F2 F3 FA 

FAlUfRr 

nmi 

701A 

0.375 

5*00 

80*0 

20 

2*00 

11.44 

10*17 

10.75 

10.6? 

11.59 

107 

111 

138 

159 

IV 

702A 

1. 000 

5*00 

45*2 

10 

1.77 

11*88 

10.58 

11.54 

10.57 

11.78 

5 

4 

4 

7 

IV 

703A 

1.000 

5.00 

78.5 

10 

1*77 

12.42 

10*76 

12.56 

11.68 

12.59 

2 

3 

5 

6 

LV 

70AA 

1,000 

5.00 

45.2 

50 

1.77 

12p59 

11.22 

11*52 

10,87 

12*12 

489 

470 

490 

491 

IV 

705A 

1.000 

5.00 

78.8 

50 

1.77 

11.75 

10*54 

11.03 

10. 3t 

11.88 

585 

588 



$ 

70f#^ 

J.OOO 

5*00 

A3*2 

10 

2*01 

11*72 

10*58 

11.07 

10.25 

11.98 

5 

4 

8 

9 

IV 

7Q7A 

1.000 

5*00 

78*8 

10 

2.01 

11.27 

10.17 

11.00 

10.10 

11.63 

2 

5 

8 

7 

IV 

70BA 

1.000 

5.00 

A5»2 

50 

2*01 

11.75 

10*17 

11.15 

10.51 

11.88 

588 

587 

5«5 

405 

IV 

70^A 

1.000 

5*00 

78.8 

50 

2*01 

11*59 

10.15 

10*79 

10.00 

i2*02 

252 

255 



5 

710A 

1.000 

5.00 

60*0 

20 

2.00 

11,41 

10.49 

11*10 

10.86 

12*01 

585 

584 



5 

7UA 

1.000 

5.00 

80.0 

20 

2*00 

10.74 

10.17 

11,02 

10.51 

11,48 

508 

510 

517 

518 

IV 

712A 

1.000 

5*00 

60.0 

20 

2*00 

12*09 

10.98 

11.91 

11.45 

12.40 

172 

175 

236 

241 

IV 

7UA 

1.000 

5.00 

20*0 

20 

2*00 

11*88 

10*85 

11.71 

11.39 

12,16 

796 

797 



5 

71AA 

1.000 

5*00 

100.0 

20 

2.00 

11*75 

10.92 

11*70 

11.59 

12,17 

45 

47 

64 

65 

IV 

715A 

1.000 

5.00 

80*0 

0 

2.00 

11.98 

11,07 

11*88 

11.51 

12.58 

1 

5 

4 

5 

LV 

7UA 

1.000 

5,00 

60*0 

40 

2.00 

11.56 

10.38 

11.01 

10,77 

11*76 

17$ 

176 

179 

182 

LV 

7 17 A 

1.000 

5.00 

80*0 

20 

1*98 

12.24 

11*13 

11.19 

10.64 

12*25 

4 

5 

17 

18 

LV 

7XfiA 

1.000 

5.00 

60.0 

20 

2.02 

13.15 

12^02 

12*25 

11 .1? 

12.58 

25 


117 


LV 

717A 

1.625 

5.00 

43.2 

10 

1*79 

12.47 

U.56 

11*88 

11*06 

12.35 

5 

4 

7 

8 

LV 

720A 

1.625 

5.00 

76.8 

10 

1*99 

12.41 

11.19 

11.29 

10.79 

12.28 

2 

5 

4 

5 

LV 

'721A 

1.625 

5.00 

45.2 

50 

1*99 

11*04 

9*96 

11*30 

10,51 

11.88 

253 

254 



5 

722A 

1.625 

5.00 

76*8 

50 

1*99 

12.58 

11*59 

11.54 

10*77 

12,00 

2 

5 

6 

7 

IV 

723A 

1.625 

5,00 

43.2 

10 

2.01 

12*95 

11*82 

11*17 

10.82 

12*07 

4 

5 

7 

8 

LV 

72AA 

1.625 

5.00 

78.8 

10 

2*01 

15.09 

11*87 

11*74 

10.40 

12.00 

2 

5 

7 


LV 

^725A 

1.625 

5.00 

45.2 

50 

2.01 

12.45 

P ^ 

11.52 

10.98 

12.70 

292 

294 

352 

353 

LV 

726A 

1.625 

5.00 

78.8 

50 

2,01 

11.58 

10.58 

n.6i 

10.98 

12.00 

2 

3 

6 

7 

LV 

727A 

1.625 

5.00 

80.0 

20 

2*00 

12*88 

11*60 

11*42 

10*79 

12*29 

5 

7 

56 

57 


728A 

1.625 

5*00 

60.0 

20 

2*00 

15.16 

12.00 

11.66 

10,49 

12.56 

6 

165 

171 


IV 

7e?A 

1.625 

5. 00 

60.0 

20 

2*00 

12*55 

11*58 

11*41 

10*79 

12*51 

7 

8 

12 

16 

IV 


-mim 




TABLE II. - LEVELS OF CHAROE RATE AND DISCHARGE 
RATE INVESTIGATED 


Power ■upply 
unit 

Charge 

rate 

_ (awpp)* 

Discharge 
rate (C2)^*^ 
(amps) 

1 

Os 373 

1.25 

2 

1. 000 

1.25 

3 

1.625 

1.25 

4 

0-375 

3.13 

r 

1.000 

3.13 

6 

1.623 

3.13 

7 

0.375 

5.00 

B 

l.poo 

5.00 

J 1 

1.625 

5.00 


(•) 


Ae used In tquatlon (1) of text. 


TABLE III. - fJIVELS OF DEPTH OF DISCHARGE, TEMPERATURE, AND END OF 
charge VOLTAGE INVESTIGATED 





T 

(“O 

ECV < 5) 
(volts) 



43.2 

10 

1.99 

H 

(d 


76.8 

10 

1.99 

■H 

tn 

43.2 

30 

1.99 

O 

u 

c 

76.8 

30 

1.99 

u 

•1 

•H 

0 

43.2 

10 

2.01 

U4 

(U 

76.8 

10 

2.01 

m 

CM 


43.2 

30 

2.01 


76.8 

30 

2.01 

u 

w 

<01 

4J 

60.0 

20 

2.00 

4J 

CS 

0 ] 

60.0 

20 

2.00 

S 

O' 

60.0 

20 

2.00 



20.0 

20 

2.00 

OD 

100.0 

20 

2.00 

H iJ 
<d G 

60.0 

0 

2.00 

0 

60.0 

40 

2.00 

< £ 

60.0 

20 

1.98 


60.0 

20 

2.02 


(a) 


As used in equation (1) of text 


TABLE IV. - ESTIMATED COEmCIENTS Of EQUATION FOR fj FAILURE CYCLE USING 
ALL THE DATA. REDUCED M(®EL USES ONLY THOSE COEFFICIENTS 
SIGNIFICANT AT THE SI LEVEL OF SICHIFICANCE 










TABU! V, - ESTIMATEI) COEFFICIENTS OF EQUATION FOR FAItURK CTCLE USINC 
ALU THE DATA. REDUCED MODEL RETAINS ONLY T»10SE COEFFICIENTS 
SIGNIFICANT AT THE 5% LEVEL OF SIGNIFICANCE 


Coo fClc lent 

Full model 

Reduced model 

Associated 

variable 

Eatlmate 
(standard error) 

Efftlmate 
(standard error) 

Bo 

2.25(0,07)A 

2,23(0.05)* 


Bl 

-.23 <0.03)* 

.•.22(0,03)* 

Xi - (CR - 1.0) /O. 625 

62 

-.14(0,03)* 

-.U(0*03)^ 

X 2 - (DR - 3.13) /I. 87 


-.21(0,03)* 

-•21(0.03)* 

X3 - (DOD - 67,2)/19,4 


.39(0.03)* 

.39(0.03)* 

0 

0 

t 

1 

X 

Bll 

-.08(0.05) 


X? 

B12 

-.05(0.04) 


X 1 X 2 

622 

-.16(0.05)* 

-.16(0,05)* 

’‘2 

B 13 

-*09(Oe03)^ 

-.09(0.03)* 

X1X3 

®23 

-.01(0.04) 


X2X3 

633 

.03(0.02) 



6W 

.06(0.04) 


Xix/, 

B2A 

.16(0.04)* 

. 16 ( 0 .(VO'^ 

X2X4 

63 .^ 

.01(0.03) 

,| 

X2X4 

64/, 

-.21(0.03)* 

-.21(0.02)* 

xg 

Number of 




diita polntfi 




(n) 

w 

127 


standard 




error of 




eHtimatc 




<S) 

o,m 

0.287 


Squared cor- 




relation 




( r 2) 

• 802 

.767 



SlgnlflcNnt at 57. Blgnlflcanca level. 




























TABLE Vm. * lARAMETEItS FOR CAMONXCAL ANALySlS OF <2 FAILURE 
CYCLE USING COEPFICIEOTS OF FULL MODEL 



JL - 1 

1-2 

i - 3 

J - 4 

Components of 
otstlonsry 
point (X,) 

-4n67 

-2.50 

0.50 

-1.17 

Ilgtnyaltie <Aj) 

0.057 

O.OOl 

-0.092 

-0.274 

Coi^oosnts of 
•Igtnvtctors 

<«tj) 






1 - 1 

-0,21 

0.18 

0.95 

0,16 


1-2 

.67 

.58 

-.04 

.46 


i • 3 

-.69 

.62 

-.30 

.19 


1 • 4 

-.17 

-.49 

-.09 

.65 


TABLE IX. • FARAMETERS FOR CANONICAL ANALYSIS OF FAILURE CYCLE 
USING COEFFICIENTS OF FULL MODEL 



j - 1 

J “ 2 

J - 3 

J • 4 

Components of 
stations ry 
point (X,) 

-2.24 

0.07 

1.04 

0.68 

Eigenvalue (Ij) 

0.032 

-0.012 

-0.091 

-0.269 

Components of 
eigenvectors 

<eij) 






1 - 1 

-0.23 

0.15 

0.95 

0.12 


1 - 2 

.78 

.45 

.07 

.42 


1 • 3 

-.56 

.71 

-.29 

.32 


1*4 

-.15 

-.52 

- .06 

.64 




TABLE X. - LEAST SQUARES ESTIMATES OF COEFFICIENTS OF EQUATIONS FOR fj CVaES 
TILL FAILURE PITTED FOR EAai MODE OF FAILURE 
SEPARATELY USING ALL THE DATA 




r 

I': 

i: 


Coefficient 

!||HPQKQnX9H| 

Short 

failures 

Associated 

variable 



^0 

2.08(0.07) 

2.45(0.10) 


n 

-.36(0.05) 

-.07(0.03) 

Xi - (CR - 1.0)/0.625 

62 

-.26(0.05) 

-.13(0.03) 

X 2 - (DR - 3. 13) /1. 87 

63 

-.25(0,04) 

-. 10 ( 0 . 02 ) 

X3 - (DOD - 67.2)/19.4 


.50(0.05) 

.003(0. 12)'*^ 

0 

0 

CM 

1 

H 

a 

hi 



y 2 

hz 


-.13(0.05) 

X 1 X 2 

622 

-.18(0.08) 


X 2 

613 


-.07(0.03) 

X 1 X 3 

623 


.08(0.03) 

X 9 X 3 

hz 


.04(0.02) 

x| 

®14 



X 1 X 4 

$24 

.09(0.05) 


X 2 X 4 

634 



X 3 X 4 

B 44 

-.16(0.04) 

-.11(0.05) 

X 2 

*4 

Number of 




data points 




(n) 

■ 

98 

29 


Standard 


1 


error of 




estimate 




(S) 

0.349 

0.0&5 


Squared cor- 




relation 




(R2) 

.775 

.924 



(a) All coefficienta algniflcant at 57. significance level except that X 4 is 
retained regardless of aignlficancc. 


I 









TAitE XI . - USAST SQUARES ESTIMATES OF COEFFICIEOTS OF EQUATIONS FOR £4 CYCLES 
TILL FAILURE FITTED FOR EACH MODE OF FAILURE SEPARATELY 
USING ALL THE DATA (127 CELLS) 


Coefficient 

Low voltage 

Short 

Associated 


failures 

failures 

variable 


Estivate 

Estivate 



(standard error! 

(standard error) 


»0 

2.24(0.05) 

2.44(0,06) 


61 

•.28(0.04) 

-.07(0.03) 

X, - (CR - 1.0)/0.625 

E 2 

•.20(0.04) 

.13(0.03) 

X 2 • (DR - 3.13)/1.87 

63 

•.23(0.03) 

-. 11 ( 0 . 02 ) 

X 3 - (DOD - 67.2)/19.4 


.44(0.04) 

. 03(0.05) <•) 

X 4 - (T - 20) /lO 

hi 



X? 

612 


-.13(0.05) 


Hz 

-.21(0,06) 


®13 


-.07(0.03) 


023 


.08(0.03) 

X 2 X 3 

033 


.04(0.02) 

y2 

*3 

01 A 

.08(0.04) 



02 A 


X 2 X 4 

03A 

B44 

-.17(0.03) 

•,12(0.03) 

T 

Number of 




date points 
<n) 

j 

97 

30 


Standard 




error of 1 

estimate 

(S) 

0,267 

0.093 


Squared cor- 1 




relation 




(R^) 

.797 

.928 



(a) All coefflciantB •Ignlficant at 57< significance level except that X 4 is 
retained regardless of significance. 




TABIE XII, - mXDWM MKELIHOOD ESTtHAIES OF COEFFICIENrS OF EQUATIONS 
FOR <2 TILL FAILURE FmSD FOR EAQI HOPE OF FAILURE 

SEPARATELY USING ALL THE DATA (127 CELLS) 


Coefficient 

Low voltegf 
felluree 

Short 

failures 

Associated 

variable 

Estimece 
^■tendard error) 

Estimate 
(standard error) 


2.27(0.06) 

2.72(0.08) 


h 

-.26(0.04) 

-.09(0.02) 

Xi - (CR - l.0)/0.625 


-.06(0.04) 

.14(0.02) 

X2 - (DR - 3. 13) /1 .87 

h 

-..25(0.03) 

-.10(0.02) 

X3 - (DOD - 67.2)/19.4 


.61(0.05) 

-.23(0.09) 

X4 - (T - 20) /lO 





012 


-.11(0.03) 

X1X2 

B22 

-.11(0.06) 


X* 

^2 

®13 


-.08(0.02) 

X1X3 

023 


.07(0.02) 

X,X3 

033 


.00(0,02) 

4 

01.', 



X1X4 

024 

.18(0.05) 


X2X4 

834 



X3X4 

644 

-.13(0.02) 

.05(0.03) 


0 

0.30(0.02) 

0,067(0.01) 



TABI£ XIII, - MAXIMUM LIKELIHOOD ESTIMATES OF C(KFFICIENTS OF EQUATIONS 
FOR £4 CYCLES TILL FAILURE FITTED FOR EACH MODE OF FAILURE 
SEPARATELY USING ALL THE DATA 
















TABU XIV. - UAST SQUARES ESTIMATES OF COEFFICIEHTS OF EQUATIONS FOR £2 CYCLES 
TILL FAILURE FITTED FOR EACH HOOE OF FAILURE SEPARATELY 
USIMC EDITED DATA CET 


CQtffici«nt 

Lon vplU|t 
failures 

Short 

failures 

Associated 

variable 

EatlMta 

fstandar<LorrprV 

Estiaate 
(atfndard error) 

^0 

2.09(0.06) 

2.46(0.09) 



-.31(0.04) 

-.09(0.03) 

Xi - (CR - 1.0)/0.625 

02 

-.20(0.04) 

.12(0,03) 

X2 - (DR - 3. 13) /I. 87 


-.23(0.03) 

-.09(0.02) 

X3 - (DOD - 67.2)/19.4 


.55(0.05) 

. 00 ( 0 . 11 ) 

X 4 - (T - 20) /lO 

hi 



x2 

*1 

012 


-.11(0.05) 

X1X2 

022 

-.16(0.07) 


X| 

013 


-.08(0.03) 

X1X3 

023 


.07(0.03) 

X9X3 

033 


.04(0.02) 

^3 

014 



X 1 X 4 

024 

.15(0.05) 


X2X4 

034 



X3X4 

«44 

-.15(0.04) 

-.11(0.04) 


Misobflr of 




data points 




<n) 

95 

28 


Standard 




error or 




eetiioate 




(S) 

0.308 

0.101 


Squared cor^ 




relation 




(r2) 

.81^ 

a900 



(a) All cpefflcientB aigniflcant at 57. significance level except that X 4 is 
retained regardless of significance. 



TABLE XV. - LEASE SQUARES ESTIMATES OF COEFFICIENTS OF EQUATIONS FOR CYCLES 
TILL FAILURE FITfED FOR EAQJ MODE OF FAILURE SEPARATELY 
US INC EDITED DATA SET 


CoefflcUnt 


Lev voltage 
failure 


Eatinato 
standard error 


2.24(0.05) 

-.24(0.03) 

-.16(0.03) 

-.21(0.03) 

.47(0.04) 


Short 

fflUure 


Estlinate 
atandard error 


2.46(0.09) 
-.09(0.03) 
.12(0.03) 
-.09(0.02) 
.00(0. U)<«) 


Aaaoclated 

variable 


- <CR - l.Q)/0.625 
X2 - (DR 3. 13) /I. $7 
X3 - (DOD - 67.2)/l9.4 
X4 • (T - 20) /lO 


-.19(0.06) 


-.11(0.05) 

-.08(0.03) 

.07(0.03) 

.04(0.02) 


Number of 
data pointe 

(n) 

Standard 
error of 
estimate 
(S) 

Squared cor- 
relation 


.13(0.04) 

-.16(0.03) 


-.11(0.04) 



(a) All coefficients significant at 57. significance level except that X 4 Is 
retained regardless of algnlflcance. 






















uhiz XVII. - mmm liwuhood estimates or coEmcxENr$ or equations 

FOR CYCLES TILL FAILURE >^1TTEP FOR EACH HOPE OF FAILURE 
SEFARATELY USING EDITED DATA SHEET 


Cocfficliiit 


Lov vpUi|t 
£«llur«n 


EitittiU 
itsndird crro 


Short 

fAilMron 


Eft IBM te 
•tandard arror 


2.37(0.04) 

2.74(0.07) 

*.21(0.03) 

-09(0.02) 

-.02(0.03) 

.14(0.02) 

-.23(0.02) 

-.10(0.02) 

.52(0.04) 

-.22(0.09) 


-.08(0.03) 

-.12(0.05) 

-.09(0.02) 


.06(0.02) 

.00(0.02) 

.16(0,04) 


-.14(0.03) 

-.05(0,03) 

0,23(0.02) 

0«066(0.01) 


Aaaociatad 

variabla 













ACTUAL DCPTH OF DISCHARGE % 

Figure 1, - Predicted cycles untU f 2 failure vs actual depth 
of discharge for five temperatures. Oata from reduced 
model coeUideiils of VL Char^ rale oT LO amp. 
discharge rate of 3. 13 amp. and end charge voltage of 
ZOO volts. 
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ACTUAL DEPTH OF DISCHARGE. % 

\ ~ *’'*'’**^ cycles until f, failure vs actual dmth 
model coefficients of TaWe VL Charoe rate of 1 n IZ 

temperature of 2CPC. endof chargeZlIageof ZOo X 


•10 


0 


90 


40 


10 20 
TOPERATURE, ®c 


h *»««« « temperature 

fc^imrw ^Khaige rates. Data from reduced mode! 


CYCLES UNTIL FAILURE 
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DISCHARGE RATE imp 

Figure E - SUrKSirdlzed residuils from predicted 
equetlon using least squares for cycles until I»r 
N ftili^re. (Data based on reduced model 
of TiMe X ) Stan^rdized residuals defined as 
• Symbols plotted Indicate number of 
points (tserved at those co-ordinates (O de- 
notes a single value, 2 through 9 indicates that 
many values, d Indicates 10 or morel 
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Figure 7. - StarHJardfzed resWoals frora predicted 
equellon for cycles until In* vollage failure. 

(Data based on reduced eijoalion of Tawe X I 
Standardized residuals defined as (yi - JiVs. 
Synbols pished Indkale number of ^Jnls ob- 
served at those co-ordinates I O denotes a sLnjf! 
value. 2 through 9 Indicates that many varl- 
abtes. O indicates 15 or morel 
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Fljur« & - Stondinllnd ratidttilt fnxn prcdlcM iqiMlIm 
using ImsI squires for cycles unlll low voKige fy fillures. 
lOili Used on coofllcitnis of TiOli XI. I Slmdordtztd re* 
sMuils <Mlno4 is |y • fVs. Symbols ploMod Indlcito num* 
bir of polnis obsoreod at llwso co^lnalos ( o dinolos 
a single value, 2 through 9 Indicates that many values, 
□ Indlales 10 or morel. 


Figure 9. • Residuals vs charge me from the 
masimum likelihood estimates for the f| shorting 
failures, (Data based on all tlw data and coeffi* 
clents of TaMt XIIL I Residuals defined as ob> 
sarvid login 114 short failure cydel minus estl* 
mated location parameter of extreme value distrl* 
butlon. Symbols plotted Indicate number of points 
observed at those co-orldlnates ( o denotes a 
single value, 2 through 9 Indicates that many 
points). 
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INVERSE EXmEME VALUE PROBABILITY SCALE 



Figurt Id • Prohiblllty plot of rvs)(iuit$ froni mixlmum tike*^ 

Mhood fit of cycles until f 4 short failure using ail the data 
plotted using extreme value prohablllty scale* Symbols 
plotted Indicate number of points obser^ at those coordh y 
natas ( O denotes a single value, 2 through 9 Indicates 
that many points)* 


Figure M ^ Probabllty plot of residuals from maximum like* 
llhood fit of cycles until short failure using all the data 
plotted using normal probability scale. SymWs plotted 
indicate number of points observed at those coordinates 
( o denotes a single value, 2 through 9 Indicates that 
many points)* 


CYCLES UNTl FAILURE 




Fiqure 11 - Expected cydes unlH fj tenure lor «ch mode se^- 
rately vs ictiHl deplh o( discharge for four teiripereluret Data 
toSd on TaMeXVL U«s char^Trtr^^^^ 

discharge rate of 1 15 amp. and end ef charge vdtage of 1 00 
voUs. 


iL 



CYCLES UNTI. LOW VOLTAGE FAILURE 
CYCLES UWTl SHORT FAILURE 





CYCUSUNTIU, FAILURE 



r^ure 15. - EqwcM cycles unUf f, miurt for nch mode 
seperalety isi function of »ctujl iMpth oTdiKfitrae lbr 
three dischatse rates. Data based on coefficients of 
We XVL Uses charge rate of 1 Damp, toipcrsture of 
2Cr C. and end of charge witage of 2, 00 wits. 





